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ABSTRACT 



A study of single-span initially bent columns is made and extended 
to include two-span columns. The axial force in these columns results 
from heating (or an equivalent) and from bowing. A snap-through 
phenomenon is considered in the analysis of the single-span and 
observations are made concerning snap-through in the two-span column. 

A digital computer program is developed to carry out the analysis and 
a description of the program enables its use by anyone familiar with 
the FORTRAN language. 
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I, INTRODUCTION 



A, GENERAL REMARKS 

The subject of initially bent columns and shallow arches is 
treated extensively in the literature. Hsu [Ref. B] studied the 
dynamic stability of shallow arches against snap- through ; Hoff and 
Bruce [Ref. 5] analyzed dynamically the snap-through of shallow 
arches; Fung and Kaplan [Ref. 2] studied the buckling of low arches 
or curved beams of small curvature; Timoshenko [Ref. 8] has also 
discussed the problem; Huddleston [Ref. 7] studied the effect of 
axial strain in elastic columns; Gjelsvik and Bodner [Ref. 3] 
approached the problem on the basis of energy criteria. There 
have been other similar studies made. However, these treatments 
of the problem have all involved only a single span. Brock [Ref. 1] 
studied the single .span using Fourier analysis in such a way as to 
provide for assembly into a multispan continuous column, and gave 
conditions for the solution of multispan cases. However, he did 
not actually arrive at numerical solution in any particular case, 
only suggesting how the computations might be carried out. There- 
fore, it was the objective of this thesis not only to study the 
single-span column in greater depth, but to develop its analysis 
in a manner such that it could most easily be extended to a multi- 
span problem, and actually to make this extension to the case of a 
two-span column. 
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B. SCOPE OF WORK 



The theory of the single-span problem is first developed for a 
column which is initially nearly straight. This development takes 
into consideration an axial strain induced in the column and also 
allows for the external application of end moments. It is assumed 
that the configuration, with and without load, lies in a single 
plane. The theory is then extended to include the two-span column. 
A digital computer program is developed to carry out the analysis, 
and a description of the program enables its use by anyone who is 
familiar with the FORTRAN language. 

Analysis beyond the two-span problem is not carried out in this 
thesis, although it is envisioned that multi-span problems can be 
investigated in much the same manner as described herein. 

C. NOTATION 

A - cross sectional area 

a - Fourier coefficients for initial configuration 
b - Fourier coefficients for final configuration 
E - Young’s modulus of elasticity 
I - moment of inertia of column 
L - length of span 
M - left end moment 
N - right end moment 
n - term index in Fourier analysis 
P - axial compressive load 
P^ - Euler’s critical buckling load 
Q - left end reaction force 
R - right end reaction force 



6 



W - lateral load in two-span column 

X - distance from left end of span 

V - initial deflection 
•^o 

y - final deflection 

$ 

e - unit strain 

- thermal unit strain 
K - curvature 

X - excess of arc length over chord length in initial 
configuration 

A - excess of arc length over chord length in final 
configuration 

7T - 3.141592 . . . 

0 - end slope of span 



7 



II. DEVELOPMENT OF THEORY 



A. SINGLE SPAN 

The study of the single span involved an initially bent column 
loaded as shown below: 

y 




Figure 1. Loading on single span column 

The end supports are frictionless hinges which do not move from 
their original positions. The moment loading shov;n is applied by 
external agencies. The initial, no-load configuration of the column 
is presumed to be given by the equation 




where here and later the summation is with respect to n and goes from 
n=l to n- ^ . For computational purposes the upper limit is taken to 
be a finite number, which v;ill be discussed in a later chapter. 

Next suppose that the column is uniformly heated so as to provide 
thermal strain. (It should be noted that an equivalent situation 
would result from moving the ends of the column tov;ard each other. 
However, throughout, we will speak in terms of heating the column 
rather than displacing its ends.) The new configuration can be 
presumed to be represented by the equation 
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The immediate purpose is to find the and thereby the new 

configuration . 

There is an axial compressive force P set up in the column such 
that the unit strain is 



e - P/AE 



(3) 



where is the thermal strain , A is the cross section area, and E 
is Young *s modulus of elasticity, (Throughout, we presume that Hooke’s 
law is satisfied.) 

The excess of arc length over chord length can be adequately 
represented for the initial deflection by the approximate relation 



o 2 



dy y 2 2 

dx = r-r- In a 



dx 



4L 



(4) 



0 



Likewise for the new configuration 



4L n 



(5) 



Thus , the change in length due to heating is 



X -X 






a 2) = Le„ - PL/AE. 
n T 



( 6 ) 



The bending behavior of the column must also be accounted for. By 
static equilibrium 

Q = -R = (N-M)/L. (7) 

The bending moment at any point is given by 



B.M, = QX + M - Py, 



( 8 ) 



and this accounts for the change in curvature 



Ak = 



2 

« d y 



d^y 



dx^ 



dx 



°= -[j) - a^) sin (2^) . 



(9) 
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According to the elementary formula. 



B.M. = El Ak . 

Thus, with a little manipulation, we find 
^ [n^P (b - a ) - Pb 



(^) 



c n n 

where we have introduced the notation 



P = . 

c 



+ Qx + M = 0, 



This is the Euler critical load which appears naturally in this 
calculation. To deal successfully with Eq. (11) it is useful to 
express Qx and M as Fourier series in the range 0 x ^ L. 



Qx = (2QL/'^) (1/n) sin 



niTx 



M = (2M/TT) PCI - (1/n) sin 



. /nnx \ 

[ir I • 



Thus , 



0 = z<(^l) 



n+1 



l^\ + Ci-(-i)"] ( — \+ n^P (b -a )-Pb \ sin ( 2 ^ 
lnir) \nTry cnn n/ \^L 



which will be satisfied if each of the terms in braces vanishes. 
Accordingly, we set 



njnc ^niry yj(n‘^p -p) 



, n- 1,2,... 



However this equation still involves the unknown, P. Equation (6 
be written in the form 

P = AEe^ -T\^AE/ni? I n^Cb^^- ) . 

Rewriting Eq. (17) once more in the form 

f(P) = 0, 

where 



( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

can 

(17) 

(18) 
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( 19 ) 



f(P) = PL/AE + (7r2/4L)5:n2 (b ^ - a ^) - Le^, 

n n T 

we seek the zero (or zeros) of the nonlinear function f(P). Because 
of nonlinearity 5 computational difficulties were encountered in find- 
ing the significant root. Of several methods attempted , the most 
suitable was the Newton-Raphson iteration scheme. A study of the 
curve f(P) revealed that the method would usually converge on the 
first root if the iteration was started at a value just slightly 

smaller than P . 

c 

The value of P thus being known, Eqs. (16) can then be solved for 

the coefficients b . 

n 

Figure 2 shows the initial and final deflection curves of a column 
computed using this analysis. The final deflection shows the column 
after having undergone thermal strain. The end moments were zero for 
this simple case. 

Further analysis of the single span problem will be given in 
Section III of the thesis. 

B. TWO SPAN 

We now consider two such spans, the column being joined at the 
intermediate support so as to maintain continuity of slope and moment 
as well as of deflection. In general, different axial loads will be 
acting in the two spans. As before, the problem is to determine the 
final configuration after heating, the initial configuration having 
been specified. 

The obvious method of analysis is, in effect, to consider two 
adjacent single spans, each of which is as shown in Fig. 1 and is 
analyzed as in the preceeding subsection. The adjoined spans are as 
shov/n in Fig. 3. Separating them and considering their mutual 
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FIGURE 2. SINGLE-SPAN COLUMN UNDERGOING THERMAL STRAIN 



interaction and the continuity at their junction, we have (see 
Fig. 4) 



Ni = 



( 20 ) 



0 = 0 . 

R1 L2 



( 21 ) 



If the common value can be determined, the problem is 

solved in the sense that each span can be now analyzed separately. As 

a matter of fact, in the solution for = M^, we will actually have 

determined P, and and the Fourier coefficients b for each span. 

12 n ^ 




Figure 3. Loading on two-span column 





Figure 4. Separation into two single-span columns 



The left and right end slopes of each span may be found by 



evaluating dy/dx at x = 0 and x = L, respectively. 

nTTb 

nTTb 



( 22 ) 

( 23 ) 
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A method of solution for the two span problem is outlined as 
follows. First assume a common value for and (This will be 

called the junction moment.) Then analyze each span separately, using 
the analysis of the preceding section, and determine the slopes at 
the junction. Compare the two slopes and find the difference between 
them, which will be called error. Next repeat the above for another 
value of moment and find the error. A plot of error versus moment 
will indicate the value of junction moment that will cause the error 
(in slope) to be zero. At this point, correct values for the Fourier 
coefficients b^ will be known for each span so that the final config- 
uration can then be computed and plotted if required. 

If the plot of slope error versus junction moment is nearly 
linear, a slope intercept method of iteration may be used to determine 
where the error goes to zero. This method simply takes two points, 
gets the equation of the line passing through these points, and then 
solves for the zero of this linear approximation. This is repeated, 
using a new point and one of the old points, until the desired 
convergence is reached. If the plot of slope error versus junction 
moment is significantly nonlinear (this occurs if the axial compressive 
load is a large fraction of the Euler load), then the slope intercept 
method is no longer satisfactory and a method such as successive 
decimation must be used. In this context, decimation means recalcu- 
lating over an interval of interest, using increments one tenth as 
large as previously. The program decimates four times so that the 
final increment is .0001 times the initial increment. This method is 
well known and will not be described here. The accuracy of results 
obtained by using this method depends upon the size of the increment 
used in stepping. The smaller the increment, the greater the accuracy. 
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The question here is where to start the iteration. If a good 
approximate figure is known, then the increments may be small. 

If not, the increments must be larger. Unfortunately there is no 
quick way of determining a good starting approximation. 

C. INITIAL DEFLECTION 

For the single span problem, any initial configuration may be 
specified subject only to the condition y(0)=y(L)=0. If the 
configuration is specified in the form of a polynomial in x, the 
axial distance from the left end, the procedure to be described 
later quickly determines as many of the Fourier coefficients a^ as 
may be desired. 

However, for the two span problem the necessity for continuity 
places extra requirements on the initial deflection. Not only must 
the deflection be zero at each end of the spans, but the slope and 
curvature must be continuous at the junction. Accordingly, such a 
two span configuration is obtained as follows. An arbitrary config- 
uration (conveniently in the form of a polynomial) is assumed which 
has zero displacement at the left of span one and at the right of span 
two. (see Fig. 5). Then, regarding this as a simply supported beam, 
sufficient lateral load is applied at the Intermediate support to 
cause the displacement to vanish at this point. The corresponding 
deflection curve is calculated and its superposition with the 
originally assumed configuration gives a result which now has zero 
deflection at each of the three support points. In general, the 
column will thus start out with non- zero bending moment except at 
the extreme ends. 
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Figure 5, Original specified configuration 



Figure 5 illustrates the originally specified configuration 
before the intermediate constraint is enforced. It is not dif- 
ficult to show that the necessary lateral load is given by 



Force = W 



Sa^+L^) EIA 

T ^ 

Li 



downward . 



(24) 



It can also easily be shown that the deflection of a beam due to a 
downward force W at point B is represented by 
WL„ 



6EI(L +L ) 

^ L 



3 , 

(L^+Lp x-x - Lj X + j— 



<x - L > ' 



where the < > symbol has the meaning 



(25) 



if X < L. 



< X - L^>- 



(x-Lj^) ^ if X 2. • 



Thus, by superposition, the new configuration is 
WL„ 



y = f(x) - 



6EI(Lj+Lp 



(L^+L ) 

3 T 2 - 12 



(Lj^+L^)^ X - x^ - L^'^x + 



<x - L >■ 



For the first span the expression for the configuration is 



WL^ X 

^ “ 6eTTl^^+lP' 



(L^+L^)^ - x^ - 



(27) 



where 0 < x '^ 12 .’ 



(26) 
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For the second span 



y = f(x5’«+L^) 



WL^Cx^+L^) 

SEKL^+L^) 






(x*+Lj^) ^ 



l } + 



... 

L^(x*+L^)^ 



( 28 ) 



where x* = distance along span two from its ovm left end. 

Now the Fourier coefficients can be found for each span. If f(x) 
is a polynomial, the program described in Appendix A makes the deter- 
mination. Thus we have a deflection satisfying the necessary initial 
conditions for the two span problem. 



o 
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III. DISCUSSION 



A. ANALYSIS OF THE PROBLEM 

In analysis of the single-span problem, the effects of the 
application of end moments were investigated. This means that 
end moments were externally applied at one or both ends of the 
column with or without the presence of axial strain. It was 
found that a snap-through phenomenon could be observed using 
this analysis. That is, by increasing one or both end moments, 
the configuration could be seen to ”jump^’ from one configuration 
to a distinctly different configuration with only infinitesimal 
variation in moment. Figure 6 illustrates this phenomenon by 
showing a column that is first heated up, and then subject to a 
left end moment. Curve 1 represents the initial configuration 
of the column. Curve 2 is the configuration of the column after 
undergoing thermal strain. The remaining curves show the config- 
uration under the action of several different values of end moment. 
It is seen that the first few 1000 inch pound increments of applied 
moment cause only slight changes in configuration, but that the 
change from 6000 to 7000 causes a relatively large change which 
has the nature of a jump, or snap-through phenomenon, following 
which addition of Increments of applied moment again result in 
only minor adjustments of the configuration. 

Figure 7 further illustrates the phenomenon by showing a plot 
of centerline deflection versus time, where time is used in the 
sense that first thermal strain is increased uniformly, with no 
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FIGURE 6. SNAP-THROUGH PHENOMENON IN SINGLE-SPAN COLUMN 
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end moment , and then thermal strain is held constant and the end 
moment is increased at a uniform rate. 

It may be noted that a similar analysis could be made for a 
column under lateral loading instead of end moment loading. 

It will be of interest to explain why the phenomenon of snap- 
through is considered to be of great importance. In a piping system 
which is heated up, considerable energy is stored in the form of 
strain energy. If snap-through or a similar phenomenon could take 
place 5 some of this energy might be suddenly released and have a 
catastrophic effect on the piping system or adjacent installations. 

In the case of a single-span straight-line system, if the ends 
should be pinned, so that no moments may be applied, it is clear 
that snap-through cannot take place. If the ends are fixed to 
heavy equipment, such as pumps or vessels, the restraining influence 
of these items would indeed result in end moments being applied to 
the pipe, but it is difficult to envision these moments actually 
resulting in snap-through. 

However, in the case of multi-span straight-line configurations, 
at the intermediate supports adjacent spans exert moments on each 
other, and it is not immediately clear that as the system is heated, 
perhaps nonuniformly , these intermediate moments might not increase 
in such a way as to cause snap-through. The intuitive feeling that 
this would not take place was insufficient basis for asserting that 
it could not happen. Accordingly, considerable attention was focused 
on this possibility. It was only after studying many different initial 
configurations and thermal loading sequences that a reasonable under- 
standing was achieved of the restraining influence which one span has 
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upon an adjacent span and of how this serves to prevent the phenomenon 
of snap-through, even if one end of a multi-span configuration is 
subjected to increasing external moment. 

Clearly, the simplest multi-span configuration is that having 
two-spans. Accordingly, this case was studied in detail and computer 
programs were written for its analysis. One of the cases which was 
studied is illustrated in Fig. 8. In this case a negative moment 
of increasing magnitude was applied at the extreme left causing the 
column to go through the succession of configurations shown. (Inci- 
dentally, although initially there was axial compression in the column 
due to initial heating, by the time the last configuration shown was 
obtained, the axial load in the left span had become tensile.) 

This example provided the insight to see why snap-through, at 
least in any sense similar to that for the single-span column, is 
not possible for two-span columns. The original, incorrect, view 
was that moment loading at the left end of a two-span column was, 
in effect, si.mply a way of applying moment to the left end of the 
right-hand span; increasing the leftmost moment sufficiently v/ould 
eventually apply a very large moment, at the junction point, to the 
left end of the right span. This is true enough, but the continuity 
between left and right spans does not permit the latter to snap- 
through as would be the case if it were the only span and the 
junction moment were directly applied at its left end. In effect, 
the left span acts as a stabilizing influence to prevent snap-through 
of the right span. 

The situation may be visualized more clearly by imagining that 
the continuity betv/een the two spans is "one-way" in the sense of 
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FIGURE 8. TWO-SPAN COLUMN WITH VARYING LEFT END MOMENT 



the mechanism shown in Fig. 9. If the two spans were joined in this 
manner , snap-through of the right span would indeed take place when 
the junction moment reached a critical value. However, with the 
true continuity of a continuous column such behavior can not take 
place . 



WITH THIS PIN REMOVED, AS 
SHOWN, CONTINUIIY IS 




RIGHT 

SPAN 



Figure 9. "One-way" continuity of two-span junction 

Attempts to create a situation in v/hich snap-through takes place 
in a two-span column have proved unsuccessful. Study of the behavior 
in these cases has added strength to the conviction that such a 
phenomenon cannot take place, even under the action of increasing 
externally applied moments at the ends or at the intermediate supports. 
However, this growing conviction cannot be understood as a convincing 
demonstration that the phenomenon cannot indeed take place, perhaps 
only for certain unfortunate combinations of physical parameters. 
Accordingly, further study of snap-through in multi-span columns is 
suggested as an appropriate topic for further investigation. 
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B. FOURIER REPRESENTATION OF THE CONFIGURATIONS 



A question arises of the ability of Fourier sine series , as 
employed herein, to represent non- zero curvature at the ends of 
spans where each sine term vanishes. Certainly, the value of 
curvature given by such a series must be zero when x=0 and x=L. 

If we have initial deflections which have non-zero curvature at 
the ends of the span (or spans), our Fourier series must be able 
to represent this adequately. Also, if we have non- zero moment 
at these points (due to the application of external moments or due 
to enforcing an intermediate constraint) this corresponds to non- 
zero curvature change, which must be adequately represented by our 
Fourier series. 

Accordingly, actual plots of curvature, as determined from 
the Fourier series , were made in several cases where the value of 
curvature could be determined by other means. Figures 10, 11, and 
12 show a comparison of curvature in a case where the initial 
configuration (in the sense of Fig. 5) had uniform negative 
curvature and viheve the intermediate constraint produced a piece- 
wise linear moment change (i.e., curvature) that was zero at each 
end and a maximum at the constraint. The final curvature is the 
superposition of these, consisting of two straight segments rising 
from -.001332 to .000666 and back down to -.001332. (In Figs. 10, 
11, and 12 these straight line segments show a very small waviness 
due to the method of plotting.) These figures also show the Fourier 
series representation of curvature, using 25, 50, and 100 terms, 
respectively. It is seen that although indeed the Foxirier repre- 
sentation of curvature does vanish at the constraint points , this 
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FIGURE 10, CURVATURE IN TV/O-SPAN COLUMN, 25 FOURIER TERMS 
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FIGURE II. CURVATURE IN TWO-SPAN COLUMN, 50 FOURIER TERMS 



0005 
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FIGURE 12. CURVATURE IN TWO-SPAN COLUMN, 100 FOURIER TERf^S 



FAIRED LINES INTERSECT LEFT AXIS 
/AT A CURVATURE VALUE CORRESPONDING 
/TO LEFT END MOMENT VALUE 
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FIGURE 13. CURVATURE IN SNAP- THROUGH OF SINGLE-SPAN COLUMN 



is a very local phenomenon and elsewhere the Fourier series provides 
a good representation of the correct values. 

Another example is shown in Fig. 13 where values of curvature 
as represented by Fourier series are plotted for different values 
of left end moment in a case of snap-through. Fairing straight 
lines through the left ends of the curves gives intersections with 
the left vertical axis which correspond to the actual moments applied 
at the left. 

Probably the most convincing verification of the adequacy of the 
Fourier representation is the relative insensitivity of significant 
results (in a particular problem) to the number of Fourier terms 
employed. In a particular case of snap-through. Table I below shows 
that snap-through occurred in the same moment interval for all values 
of n down to as lov; as 3. It is interesting to note in Table II that, 
even for n=l, the axial load obtained is in very good agreement with 
the load obtained using 100 terms. 



Table I. Values of end moment between which snap-through 
occurred, for different values of n. 



Moment (in .lb.) 


n=100 


o 

LO 

II 

C 


LO 

CM 

II 


II 

O 


II 

cn 


CO 

II 

C 


n=2 


1 — 1 
II 

c 


60450 

60460 

60470 


9 


• 


• 


• 


* 


€ 


• 





Note: Dots indicate the moment interval in which snap-through 

occurred. The tabulation is for a particular problem 
described on page 31. 
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Table II. Value of axial load in pounds for different 
values of n and M. ■ 



Moment ( in . lb . ) 


60450 


60460 


60470 


n=100 


38489.000514 


38491.906680 


38489.899526 


n=50 


38489.000515 


38491.906680 


38489.899526 


n=25 


38489.000517 


38491.906680 


38489.899527 


o 

II 

P 


38489.000550 


38491.906674 


38489.899548 


n=5 


38489.000764 


38491.906634 


38489.899684 


n=3 


38489.009402 


38491.905053 


38489.905163 


n=2 


38489.011286 


38491.009623 


38489.906358 


n=l 


38489.011298 


38491.009624 


38489.906349 



The tabulation is for a particular problem described below. 

Due to computational difficulties in the vicinity of snap- 
through, values of moment taken at smaller than 10 inch-pound in- 
crements provided somewhat erratic results. 

The results tabulated in Tables I and II are for a steel pipe 
whose cross sectional area is 2.945 in^. The length is 200 inches, 
the moment of inertia of the cross section is 5.2 in^, and Young’s 
modulus of elasticity is 30 x lO^psi. Moment is applied at only 
one end of the span. The initial configuration was similar to that 
shown in Fig. 6. 
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IV. ANALYSIS OF A TYPICAL TWO-SPAN PROBLEM 



Using the program described in Appendix A, configurations 

resulting from the application of many different combinations of 

end moments and strain on a two-span column can be obtained. The 

following cases illustrate several of the ways in which the program 

may be utilized. These examples involve a two-span steel pipe with 

2 

cross sectional area of 2.945 in . The left span is 100 inches 
long and the right span 300 inches. The moment of inertia is 5.2 in^ 
and Young *s modulus of elasticity is 30 x 10^ psi. It should be 
noted, however, that a column of any constant cross section or 
material could have been used. 

Case I - As illustrated in Fig. 14, the pipe is subjected to 
thermal strain only. Curve 1 shows the initial configuration. 
Subsequent numbered curves represent the configuration correspond- 
ing to increasing values of thermal strain. 

Case II - Here the pipe is subjected only to an increasing 
moment. Figure 15 shows the sequence of resulting configurations. 

Case III - This example shows the resulting configuration of 
the pipe (see Fig. 16) with simultaneously increasing thermal strain 
and left end moment. 
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FIGURE 14. TWO-SPAN COLUMN WITH INCREASING THERMAL STRAIN 




o 

o 
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FIGURE 15. TWO-SPAN COLUMN WiTH INCREASING LEFT END MOMENT 



INITIAL CONFIGURATION 



o 

o 
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FIGURE 16. TWO-SPAN COLUMN WITH SIMULTANEOUSLY INCREASING 
THERMAL STRAIN AND LEFT END MOMENT 
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V. RECOMMENDATIONS 



If future studies are made using the analysis set forth in this 
thesis, the follovring areas are recommended. 

The problem should be extended to include multi-span columns. 

In the analysis of two or more spans it would be convenient to have 
a better way of obtaining the junction moment or moments* Perhaps 
a method can be developed that does not require an initial approxi- 
mation. 

Another area of study would be to develop a proof that snap- 
through cannot take place or develop a clear delineation of the 
circumstances under which it could take place. The study probably 
should focus initially on t\'7o-span columns. 
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APPENDIX A - DESCRIPTION OF COMPUTER PROGRAM 

A. GENERAL REMARKS 

The digital computer program used to carry out the analysis set 
forth in this thesis is a double precision FORTRAN language program. 

It is self contained except for one outside subroutine used to draw 
the configurations. This subroutine or a similar one should be 
available in any computer library. Various write statements are 
interspersed throughout the program to allow the user to follov; the 
progress of the computations ^ and to help trace any errors. These 
statements may of course be left out of the program to lessen the 
amount of unneeded output. 

B. BASIC FUNCTIONS PERFORMED BY THE PROGRAM 

The basic program for the single-span plots the initial deflection, 
determines the axial load on the column, and plots the final deflection. 
The axial load is determined using a Newton Raphson technique on Eq. 
(19). Iteration continues until convergence is reached. 

For the two-span problem, a process is employed of repeatedly 
selecting values of the junction moment. Then .the two spans are 
analyzed and the two junction slopes are computed and compared. Then 
the correct value of the junction moment is found using the method of 
successive decimation. Although the slope intercept method is more 
accurate when it works, its unreliability causes it to be less 
desirable for use in the general program. After the junction moment 
is found, the final deflection is then plotted. 
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C. INITIAL DEFLECTION 



The Fourier coefficients for the initial deflection of the column 
are computed for the analysis previously described. The program 
accepts coefficients describing a polynomial curve which passes 
through the left support of span one and the right support of span 
two. Then the modification of this curve is computed which enforces 
the intermediate constraint; this modified curve is also in the form 
of a polynomial. This yields two polynomial curves, one for each 
span; cf. Eqs. (27) and (28). Then the Fourier coefficients a^ are 
evaluated for each span. 

1. Subroutine SINCOE 

This subroutine evaluates the integral 

7T 



GA(m,n) = y^ sin(ny)dy 



(A.l) 



which is used in the evaluation of the Fourier coefficients. The 

method involves evaluating integral number 6, page 117, of Ref. ^ 

See program listing in Appendix C. 

2. Evaluation of Fourier Coefficients 

The Fourier coefficients are evaluated according to the 

following formula: 

L 



a 

n 




sin 




n 



1 , 2 ,..., 



(A. 2) 



where L is the length of the span and n is the number of coefficients 
desired. For the left span, y, as given by Eq. (27), is substituted 
for the value of gCx), and for the right span Eq. (28) is used. 

The method carried out in the program can best be illustrated 
through the use of an example. Suppose the initial configuration f(x) 
is given by the equation 

f(x) = Ax^ + Bx. (A-3) 
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Then Eqs . (27) and (28) become 



WL, 



y = Ax^ + Bx - 



SEia^+L^) 



X 



CL^+L^)^ - x 2 - L^2 



] 



(A. 4) 



and 



“‘■2 

y = A(x.HL^) t B(x.VL^)- 



(L +L f- (x*+L 



(yL^) ^,3 

L2(x*+Lj^) 



(A. 5) 



respectively. Equation (A. 4) can be rewritten as 



y = Ax^ + Bx - Cx + Dx' 



(A. 6) 



and Eq. (A. 5) can be rewritten as 



y = A(x*+Lj^) + B(x*+L^) - C(x*+Lj^) + D(x=’-+L^)^ - Fx* 
where 



...3 



C = 



WL2(L^+2L^L2) 

BEKLj^+L^) 



(A. 7) 



(A. 8) 



and 



D = 



F = 



WL^ 

BEKLj^+L^) ’ 



W 



6EI 



(A. 9) 



(A. 10) 



For the first span, Eq. (A. 2) becomes 

L, 






fl fx"sinfi4i)dx t B fx sinf2^] 

■I \ b/ I \h) 



dx 



- C 






0 

dx + D j x3 sin 1 dx 






) 



The integrals are all of the fonn 

L, 



/ m . / nirx \ 

X 



dx 



(A. 11) 



(A. 12) 
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Making the substitution y = j this integral may be put in the form 




n + 1 „ f ^ 

P GA(m,n), 



where 




(A. 13) 



(A. 14) 



and the definition, (A.l), is used. Thus Eq. (A. 11) can be rewritten 
in the form 



a = — AP^GA(2,n) + BP^GA(l,n) - CP^GA(l,n) + Dp“GA(3,n) 

n 1 ^ 

= — ^ A. P^'^^ GA(i,n) (A. 15) 

^1 ^ 

Note, from Eqs. (27) and (28), that the maximum degree of the 
adjusted polynomials is Max(3,m'0 3 where m-’» is the degree of f(x). The 
number n of Fourier coefficients to be employed, and the integer m- 
are among the inputs to the program. Internally, m* is increased, if 
necessary, to a minimum value of 3. Then SINCOE is used to generate 
the numerical constants GA(i,j), i = 1^2, m-; j = 1,2, n. 

Thus, the coefficients are available for the evaluation of the Fourier 
coefficients as represented by such equations as (A, 15), the coef- 
ficients A, B, C, . . . , in which are obtained in. an obvious manner. 

The Fourier coefficients for the second span are found in 
basically the same manner, except that now the computations are more 
numerous due to the necessity of expanding powers of (x^HL). It can 
be seen that for the second span Eq. (A. 2) becomes quite lengthy since 
g(x) will contain many terms. Looking only at the first term in the 
example, we have, after expansion 
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(A. 16) 



< A 



^2 , 

I X‘-^ sin(^^^-^\dx + L I x* sin f 

^ ^ 2 / 1 Jo V S / 



dx 



+ • • • 



This can also be written in the form 



n L. 



<A p 






p2GA(l,n) + L 



^2pGA(0,n)|+ * 



where now 



P = 



(A. 17) 



D. VARIATIONS IN THE PROGRAM 



Different combinations of strain and end moments to be applied to 
the column can be incorporated into the program. This can be accom- 
plished easily through the use of the input cards described below. 

The program can also be altered through an additional D0~loop to shov/ 
a "time history" of the column. The strain and/or end moments can be 
incremented in the DO-loop so that intermediate configurations may be 
plotted. This is useful in looking at snap-through. The program 
listed in Appendix C has this variation incorporated. 



E. USE OF THE. PROGRAM 

This program may be used by inserting the follov/ing 21 input cards 
in the designated places at the beginning of the program. Each card 
contains only one entry. Note that these cards are made a part of the 
program deck and are not read in as data following the program. See 
program listing for example. All cards must be punched beginning in 
column seven. 
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DIMENSION GA(MPOL 5 NFOR) - integer numbers must be 

inserted for MPOL and NFOR 

MPOL = maximum degree of initial deflection function 

NFOR = number of Fourier coefficients desired 

NDEC = number of decimations in determining junction 
moment 

C 6 = initial function coefficient 
C5 = initial function coefficient 
C4 = initial function coefficient 
C3 = initial function coefficient 
C2 = initial function coefficient 
Cl = initial function coefficient 
LI = length of left span 
L2 = length of right span 
E = Young’s modulus of elasticity 

2 

W = cross section area of the column in inches 

Z = moment of inertia of the column in inches ^ 

EPSl = applied strain in left span in inches/inch 

EPS2 = applied strain in right span in in/in 

SS = applied left end moment in pound inches 

UU = applied right end moment in pound. inches 

XMOMT = starting moment for successive decimation 
iteration 

DELTA = initial increment for successive decimation 
iteration 



Notes : 

1. All inputs are floating point numbers (D format) except MPOL 
and NFOR 5 which are integer numbers. 

2. There must be an input card for each input listed, even if it 
is zero. 
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3. Each initial function coefficient must correspond to the 
power of X it is associated with. That is 

f(x) = C6x^ + C5x^ + C4x‘* + C3x^ + C2x^ + Clx 
(A polynomial of degree greater than six can be accommodated only by 
making internal changes in the program: these should be rather 

obvious . ) 

4. The program is dimensioned for a maximum number of 50 Fourier 
coefficients and a maximum combined span length of 400 inches. (This 
limitation is related to the fact that deflections are computed at 
two-inch intervals and are stored prior to output.) If the user wishes 
to extend these limits, the dimension statements must be altered 
accordingly. 

5. Likewise the call statement for the DRAV/ subroutine may have 
to be altered to accommodate the particular problem in question. (For 
example, to provide a desired scaling.) 

6. Where to start the iteration for the method of successive 
decimation, i.e., selection of the input XMOMT , poses a problem. 
Experience and perhaps a trial run will best produce a suitable start- 
ing point. 
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APPENDIX B - FLOW CHART 
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A AAA 



WWW 



T 

ii 

COMPUTE AND PLOT 
INITIAL DEFLECTION 






i 






y 

COMPUTE LOAD 

COMPUTE FOURIER COEFFICIENTS, b 
COMPUTE END SLOPES 







COMPUTE AND PLOT 
FINAL DEFLECTION 

1 

COMPUTE END SLOPE ERROR 
USE METHOD OF SUCCESSIVE 
DECIMATION 



y 



BCOUNT = BCOUNT+ I .0 




= 5 



( STOP ) 
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